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The energies of a pair of strongly-interacting subsystems with arbitrary noninteger charges are 
examined from closed and open system perspectives. An ensemble representation of the charge 
dependence is derived, valid at all interaction strengths. Transforming from resonance-state ionicity 
to ensemble charge dependence imposes physical constraints on the occupation numbers in the 
strong-interaction limit. For open systems, the chemical potential is evaluated using microscopic 
and thermodynamic models, leading to a novel correlation between ground-state charge and an 
electronic temperature. 
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Associating fractional charges with individual atoms 
is a habitual part of our every-day thinking about con- 
densed and molecular matter. Indeed, characterizing 
the energetics of such systems in terms of dynamically- 
evolving charges is now recognized as key to understand- 
ing the atomic-scale behavior of complex processes rang- 
ing from alloying to motor protein function to molecular 
logic gate operation Compounding the difficulty of 
the problem is that the instantaneous redistribution of 
charge occurs under the influence of strong interactions 
among subsystems of atoms and molecules. 

Historically, it was assumed that quadratic electro- 
static interactions constitute a reasonable representation 
of the energy dependence on fractional charge 0, in- 
dependent of the strength of interaction. However, the 
inadequacy of the historical assumption had been high- 
lighted by the work of Perdew et al. (hereafter referred to 
MS PPLB) 0, where a linear dependence was found for 
weakly interacting subsystems. This was found by con- 
sidering an open subsystem that was allowed to weakly 
interact and exchange electrons with a reservoir of elec- 
trons For a diatomic molecule, the weak interaction 
restriction implies that the theory is valid only at large 
internuclear separations R. PPLB considered an atom A 
as a subsystem in its neutral state with M electrons and 
energy Em that becomes anionic by fractional charge q, 
where the anionic state with M + 1 electrons has energy 
Em+i- The energy E^(q) was shown to be the ensemble 
average 



E A (q) = E M + to (E M+ i - E M ) , 



(1) 



with lo ~ wpplb = — q > 0. This result established a 
seminal extension of density functional theory (DFT) to 
fractional numbers of electrons The energy E^(q) 
is manifestly linear in q, in contrast to the historically- 



assumed quadratic dependence. A direct consequence of 
Eq. (fTJ) is that the associated microscopic chemical po- 
tential fi = —dE\(q)ldq exhibits a discontinuity with 
respect to charge q [3(. The veracity of both limits 
has been confirmed numerically by Cioslowski and Ste- 
fanov [f|. In this Letter, we construct a novel, ana- 
lytic model explicitly linking the weak (linear) and strong 
(quadratic) interaction limits, and characterize the cor- 
responding behavior of fi with respect to q. As discussed 
below, this requires the definition of interacting "atom- 
in-molecule" (AIM) subsystems and associated fractional 
charges. While these definitions are not unique 0], 
and many physically-reasonable alternatives are possible 
d, [E EE EH j our results are independent of the details 
of a specific AIM approach. 

Previous attempts to extend PPLB's formal results to 
moderate and strong interactions have recast the prob- 



lem in terms of charge resonances [12|, [13|, [14|, [15[ of a 
closed system composed of A+7Z, where 1Z represents 
the electron reservoir. In the charge resonance view, one 
supposes that special wavefunctions can be constructed 
such that the charges on the subsystems are integers 
3 EE EE EE, HE]- Interpreting PPLB in the language 
of resonances, the state designated as M for A corre- 
sponds to both A and 1Z being neutral (covalent reso- 
nance), while M + 1 corresponds to A being anionic and 
1Z being cationic (ionic resonance) , leaving the total sys- 
tem neutral. 

Extending this picture to the intermediate interaction 



regime [12, |2l| , 1Z ceases to be a structureless reservoir 
of electrons, but must instead be viewed as a structured 
subsystem B at a finite separation R from A. This sub- 
system may be an atom, molecule, or bulk material. Al- 
though several two-state valence bond models have been 
developed in this limit [H EE El El, previous attempts 
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to extend the valence bond approach to the strong inter- 
action limit have proven unsuccessful. 

A complete generalization of PPLB to the strong in- 
teraction limit can be attained by transforming from the 
resonance-state basis to an ensemble or spectral repre- 
sentation 0]. Let H denote the hamiltonian for the 
closed system AB with eigenstates {^k} and eigenen- 
ergies {E^}. Let vff be any arbitrary trial wavefunction 
with = 1. The variational energy can then be 

expressed in a spectral representation, 



(2) 



where the uJk = \{^\^k)\ 2 are occupation numbers. Now 
we express 'F as a linear combination of resonance states 
{ipi} with coefficients {c 2 }, such that "J = JZi^V'i- The 
eigenvectors Ck are the specific values of coefficients cor- 
responding to the eigenstates. The basis rotation from 
the {tpi} to the {^k} in terms of the Ck lead to occupa- 
tion numbers 



! = \ ^2 c *Ctk\ 



(3) 



These occupation numbers give rise to piecewise linearity 
in the energy when AB dissociates [H, HH . Specializing 
to a two-state model, we define the ionicity 7 as the ratio 
Ci/co, where Co and ci are the coefficients for the cova- 
lent and ionic resonances, respectively. Let 7 gs and 7 XS 
denote the eigenionicities corresponding to the ground- 
and excited-state values of Ci/cq, at some chosen separa- 
tion R 12j. The dependence of the energy E(-f) on an 
externally-imposed ionicity 7 is then given by 

E( T ,R) = E gs (R) 

+ u( r ,1 gs (R),7xs(R)) (E XS (R) - E gs (R)) , (4) 



with 



w(7; 7gs,7x S ) = 



(7-7g S ) 2 



(7-7gs) 2 + (^ 2s r) (7-7* s ) 



(5) 



This occupation number governs how the energy of AB 
changes when its ionicity is forced to deviate from its 
ground-state value (Fig. Q]). Importantly, this result is 
independent of the strength of the interaction. As is nec- 
essary physically, u> = when 7 = 7 gs , and u> = 1 when 
7 = 7 XS . Moreover, both extrema are represented, as the 
derivative of to with respect to 7 is zero when evaluated 
at either 7 gs or 7 XS (Fig. [2^). 

Equation ([5]) is expressed in terms of the ionicity 7. In 
order to complete the model it is necessary to eliminate 
this quantum mechanical parameter in favor of the phys- 
ical charge q [23| . This can be accomplished through an 
AIM decomposition 0, [ljj . The relation between 7 and 
q in the strong interaction regime is illustrated in Fig.[5}3. 
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FIG. 1: (Color) Energy of a heteronuclear diatomic molecule 
as a function of bond length and charge as modeled by Eq. Q . 
The charge dependence is linear as the molecule dissociates 
Q|, while being quadratic near the ground-state charge when 
R is near equilibrium values 0, E|. The constituent atoms 
interact strongly near equilibrium. 
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FIG. 2: (a) (Color online) Intermediate limit of Eq. ([5j. The 
extrema correspond to 7 = 7 gs and 7 = 7 XS = —1/7^. (b) Re- 
lationship between charge and ionicity in the strong interac- 
tion regime for the two-state model and the atom-in-molecule 
model of Ref . [l2| • S c i and SN^ are AIM parameters entering 
the 7 — q relation in the strong interaction regime. The gray 
zone corresponds to the physical range of q between and 1. 



It reduces to 7 = ±1 
action regime [HI, [H 



'(1 — q) in the intermediate inter- 
yielding u(q) ~_Wint(9,g gS (i?)) 

Although 



= (^(l-q gs (R))q ± ^ qgs (R)(l-q)) 2 \U 
details of the 7 — q relation will be affected by the spe- 
cific choice of AIM charge definition, we can immediately 
identify two quite general consequences of the transfor- 
mation to q. First, in specifying the 7 — 9 relation, choos- 
ing a root, and restricting the charge to the physical 
range between and 1 , some values of the ionicity are ex- 
cluded. This outcome was anticipated by Pan et al. [2~H ]. 
As a result, to(q) — in contrast to 0^(7) — may no longer be 
a proper occupation number, since it does not span the 
full range from to 1. Second, the presence of the two 
roots results in two branches in the E(q) surface (Fig.[3K; 
only one was shown in Fig. [1]), provided that the gap in 
Eq. (TJ| is nonzero. The branches are degenerate for all q 
in the weak limit, degenerate only at the integer charges 
in the intermediate limit, and completely nondegencratc 
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in the strong limit. 

As the driving force for such charge transfers, the con- 
cept of chemical potential presupposes an ability to spec- 
ify subsystems. However, this concept requires the iden- 
tification of the subsystem energies, as well as charges. 
As with the AIM charge, the subsystem energy is not 
uniquely defined. Regardless, certain general features 
are valid for any definition of subsystem energy, pro- 
vided that the total energy of the parent closed system 
is preserved. That is, for closed system energy £ab> the 
open subsystem energies £? A and must sum to Eab- 
E A + Eq = Eab- (Asterisks indicate open-system sta- 
tus with respect to energy and electron transfer.) This 
proviso is necessary because we require the stationary 
properties of the cigenstatcs of the closed system. 

Now, if -Eab is either the ground or an excited state 
energy, it must follow that the subsystem energy defini- 
tion applies to these as well, thereby defining E A and 
EXxs- Thus, for subsystem A*, its energy E* A with an 
arbitrary charge q becomes 



and thermodynamic definitions 14, [26|, 27 1. The ther 



El{q) 



EL 



(6) 



where w(q) comes from Eq. ([5]) and the exact 7 — q rela- 
tion. An analogous expression holds for B*. It is trans- 
parent that the sum of the two AIM energies recovers the 
total energy of the closed system, independent of the de- 
tails of the chosen subsystem energy definition. Likewise, 
if -Eab represents a small change from a ground-state en- 
ergy, then one concludes that —dE^/dq = dE^/dq = fx*, 
where fi* is the microscopic chemical potential. Differen- 
tiating Eq. ©, and using the simpler intermediate regime 
with the negative root of the 7 — q relation, yields 

f(q) = -(e: s -e; s ) 



x (l-2 ggs ±(l-2 g )J (1 (l ^^ ), (7) 



modynamic definition is relevant because the number of 
states in the spectral representation, Eq. (|2|), becomes ex- 
ponentially large with system size. The grand canonical 
ensemble is most appropriate here 0, 12a ]. For an en- 
semble characterized by chemical potential [i and inverse 
temperature (3 = l/k^T, the grand partition function 
is J^k ex P(l3nN£ - f3El), where N% denotes the number 
of electrons in state k with energy Et [2§|. Considering 
again a two-state model, the average electron number is 



N. 



(N 

solution for 



(^ s -^ g *s)exp(/W x * s - 
- N; s ) - (3(E* S - 

li(q) with q = N* - (N 



gs 

gs ))/(l + exp((3/i(N* s 



N* ) 

gsJ 

- E; 



E* s - E*+k B T\n( 



-If* 



N* - N* 



) whose 
O is 

(8) 



At T = K, n = (E* s - E* gs )/(N* S - JV g * s ) for all q. By 
comparing Eqs. (O and (J8j> , we observe that the T = 
K result corresponds to q gB = (Fig. [3b). On the other 
hand, a temperature of 34800 K corresponds to q gs = 1/2, 
if energy E* s -E* s = 3 cV and N* s -N* s = — 1 are chosen. 

Because there arc no other degrees of freedom with 
which to equilibrate, it is clear that T corresponds to an 
electronic temperature. In fact, this is the grand canon- 
ical (open system) an alog of the electronic temperature 
introduced by Kohn [301 ] in the context of calculating 
excitation energies using closed system ensemble DFT 
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32|. The corresponding canonical ensemble is de- 
fined in terms of a spectral representation, as here, and 
is characterized by a temperature 0. 9 is defined implic- 
itly through a self-consistent relation equating the total 
entropy of the system to the integrated local entropy of 
the system's electron density distribution, in turn defined 
through a set of state- and temperature-dependent Kohn- 
Sham equations [13, EH . There is thus a direct analogy 
between our electronic temperature for an open system 



where the R dependence and atom subscript have been 
suppressed for clarity (Fig. [3b). fi* is not defined at 
integer charges, thereby preventing the subsystem from 
transferring a full charge. This behavior comes from the 
degeneracy in the branches of E(q) at integer q and is in- 
timately related to the derivative discontinuity found by 
PPLB. At q = q gs with the negative root, fi* = 0, in ac- 
cordance with our zero of energy. Nonetheless, chemical 
potential equalization [3] does hold. 

By contrast, for strong interactions with the exact 7—5 
relationship, /i* is continuous at and 1. Smoothing of 
the chemical potential arises directly from the movement 
of the derivative discontinuities to noninteger q outside 
of the physical range of [0, 1] (Fig. [3b)- The subsystem 
can now transfer a full charge, but the driving force may 
be very large." 

To complete the discussion of chemical potential and 
in analogy to PPLB, we compare microscopic (Eq. ([7])) 
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FIG. 3: (a) cu(q), corresponding to branches of the energy 
surfaces, Eq. @, in the strong and intermediate limits for 
each root of the 7 — q relationship: iS c i = 0.25 and 5N*i = 
—0.20 (strong); S c i = and 5N* { = (intermediate). q gs = 
0.25 in both cases. Gaps form at integer charges, (b) Two- 
state microscopic and thermodynamic models of the chemical 
potential. The microscopic model is plotted for q ss — and 
1/2, and the thermodynamic model for T = and 34800 K. 
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with charge transfer, and that of Kohn et ai, even though 
they restricted their attention to polarization excitations 
of a closed system with fixed total charge. In both treat- 
ments, the temperature corresponds to deviation of the 
charge distribution away from the ground state distribu- 
tion. From a thermodynamic point of view, the excita- 
tions lead to an increase in entropy. The correspondence 
observed between T and q gs stems from the fact that 
both quantities encode information about the energetics 
of the microscopic states of the interacting subsystems. 
Equivalcntly, both T and g gs reflect the strength of the 
coupling between a subsystem and a reservoir. 

In conclusion, we have shown how charge dependen- 
cies in a pair of subsystems can be described by appeal- 
ing to an ensemble representation of the closed system 
energy, even in the strong interaction limit. The an- 
alytical form of the charge-dependent system energy is 
determined through a basis rotation involving the reso- 
nance states that embody charge transfer between sub- 
systems. When a decomposition method is applied to 
closed system eigenenergies, charge-dependent, open sys- 
tem energies can be defined, leading to a microscopic 
model of the chemical potential. When the subsystems 
interact strongly, the derivative discontinuities in the en- 
ergy found by PPLB move to noninteger charges that are 
outside of our physically allowed range, leading to con- 
tinuous chemical potentials at the integer charges. When 
compared with a thermodynamic model of the chemical 
potential, a correlation is observed between the ground- 
state charge and an electronic temperature, thus defining 
complementary measures of the interaction strength. 

The ensemble variational energy defined here may be 
regarded as the wavefunction predecessor of the excited- 
state density functionals of Gross et al. 01 . However, 
the formal connection of our strongly-interacting, open 
system results to density functional theory remains an 
open question, requiring the construction of a variational 
principle for open system excited states. To the best of 
our knowledge such a variational principle has not yet 
been established. 
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